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Abstract 

In this paper, the fundamentals of the structural 
diagram method [1] for Distributed Parameter Systems 
(DPS) are presented and reviewed. An example is given to 
illustrate the application of this method for control design. 
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I. Introduction 

For lumped systems, practical control methods are 
developed based either on the transfer function or the state 
space representation of the system [3], Although there are 
many inherent connections between these two methods, there 
exists many differences between them in terms of 
mathematical analysis tools, design methods, and 
engineering interpretation, etc. For DPS, however, control 
methods are predominantly based on the state space 
representation of the system [2). It seems that this approach 
parallels that developed for lumped systems, but it is 
difficult to understand for practicing engineers, and to 
implement a physical controller without an intensive 
functional analysis background [2]. The structural diagram 
method for the DPS (the counterpart of transfer function 
method for lumped systems), on the other hand, should be 
easier for practicing control engineers due to its simplicity, 
convenience, and clarity of representation of the systems. 
Also, the popularity of this method should be enhanced by 
the fact that control engineers are already familiar with 
complex variable methods from their experience in dealing 
with lumped system. 

In this paper, the fundamental concepts of the structural 
theories arc reviewed in an attempt to explore the 
methodological value of a generalization of the structural 
theory of DPS for practical control design. 


transform of this signal w'ith respect to time is w(x,s). 
2. Distributed block 


U(x,s) WW(y,x,sl 


*Y(y,s) 


In order to describe the input/output relationship of a DPS, 

a distributed block, an analog to the transfer function of a 

lumped system, has the form HYyxs), x£Dj t yeD 2 for a 

stationary process (the coefficients in a PDE do not change 

with time), where Dj and Di are the domain for the distributed 

input and output signal respectively. Therefore, the 

relationship between any two places in die space is reflected 

by inserting in their position in HYy.x.J). This is then the 

Laplace transfer of the Green’s function of the given DPS. The 

relationship of the input and output is not simply a product of 

the input signal and transfer function as in the lumped system, 

but instead, a composition of the two, which is defined by: 

r 


Y(y,s) = \V( yjc,s)SU(x,s) = 


W(yjc t s)U(x t s)dx 


Jd 


2. Connecting blocks 

A complicated DPS may involve more than one block. 
The operations among them then need to be clarified. The 
following three kinds of operations often arise in control 
design. 

Connections in series. 


U(x.s) — > |>Vl(.xI,.\7] Kv2;y,xl,sl — £> Y< y.s) 

If two blocks are connected in series as shown above, the 
overall block is equivalent to: 

\V( yj,s) = \Vl(xijc t 5)&\V2(xi,x,s) = \V](xu t s)\V2(yjj t s)dxi 

JD 


This is again, different from the lumped system where the 
overall transfer function for two systems connected in series 
is simply the product of two individual transfer functions. In 
addition, the operation of composition is non-commutative in 
general. If both block functions are known, the overall block 
rfvis obtained after conducting the conventional integral 
# operation. 

Closed loop system. 


II. Structural Diagram Method 

The structural diagram method developed in [1] is 
briefly reviewed here. 

1. Distributed signal 

In DPS, a distributed signal depends not only on time 
but also on certain spatial variables. It then has the form 
xeD y where D is a spatial domain. The Laplace 


U(x,s) 


TJ 


|Wl(y,x,s) 


-J\V2(x,y,s) 


-►Y(y,s) 


Based on the block operation introduced before, it is easy to 
show that the closed loop transfer function W(y,x,s) has the 
following general form: 



W(yj.s) = W 2 ,(yjc,s)(S>W(yjc.s) + Wi(yj.s) 

W 2 l(yJ.s) = W I (ypc.s)®W 2 (yJ.s) 

For special types of Wl(yjcj) and W2(x,y.s), W(y.x,s) may 
be simplified greatly. 

Lum ^ C p 0 Howing figure shows that the interaction between a 
DPS and a lumped parameter regulator is represented by two 
§ functions. The closed loop transfer function from U(x,s) to 
Yfa.s) is simplified to: 

WHajc.s) 

W(ajcj) 


l-H(s) 


Wl(ajcj)dx 


JD 


U(x,s) 



4. Transfer function of DPS 

The transfer function of DPS is closely related to the 
attached boundary condition, as weU as given PDE. This is 
what differentiates a DPS and lumped system. The transfer 
functions of many PDEs with popular boundary conditions 
are given in [1], For a general PDE with arbitrary boundary 
conditions, one needs to find the new transfer function by 
solving the given PDE with its boundary condiuons. 


The system written in standard form is. 

dO(x.t) jdQ(xfL +Q(x.t) = u(t)+S(x)d(t) 
Bt dx 


d(s) 



The open loop transfer function W(ys,s) for the given plant 
and boundary conditions is easily derived or found by looking 

in the table in [1] to be 

W(yx.s) = 

The closed loop transfer function from d(t) to Q(1 4) is 

HW— 


rl 


1 +H(s) 


\V(ljc.s)dx 


(s+l)e 




( s+l)+II(s)0-e <”’>) 

The dosed loop transfer function from Q*(t) to Q(l.t) is 


H(s) 


W 2 (s) = - 


\V(1 jc.s)dx 


1+1 1(s) 


W(]jc.s)dx 


' (s+U+lHsXJ-e^) 


III. An Example 

The following figure shows a furnace of length 1=1 
which heats a continuous snip of material passing m.oug i 
the furnace with speed v=i. 


Furnace 


jFucl u(t) 

\P7ZZZZ] ▼ r ZZZZZ] .Material 

YA t-d being healed 


Let the temperature field of the furnace be uniformly 
distributed in space. Suppose the temperature distribution in 
the heated strip isQ(x.t). Then for a normalized system, u(t) 
and Q(x,t) are connected by a first order linear partial 

differential equation 

BQQcj) JWxdl +Qfxj) = u(t). 0<x<l, t>0 
Bt dx 

w-ith boundary condition Q(0,t)=d(t). 

Two specific control tasks can now be considered. One 
is to regulate the temperature at the output of the furnace. 
This is considered first below. The other task is to heal the 
material to a specific spatial temperature distribution. 

For the first control task, the control equation has the 
form U(s)=W(s)[Q(1.0-Q*(l't)]- Then the control diagram is: 


The response of Q(1.0 * s l * lcn 

0(1, s) = \V,(s)®d(s)+W 2 (s)®Q*(l ,s) 

The characteristic equation cf the system is 

A(s) = (s+l)+H(s)(l-e^ 1 >) 

based on which, we can carry on control analysis and design. 
All this can be done using conventional control techniques, 

such as root locus. Bode plot, etc [3]. 

Consider a proportional controller where H(s)=k. ine 
poles of the closed loop system arc shown in Fig. 1. Notice 
that there is an infinite number of open loop zeros and one 
open loop pole at s=-l that is cancelled by the zero at the 
same place. It is seen that for all positive k, the closed loop 
system is stable. The closed-loop poles for a general PID 
controller are shown in Fig.2 where 
H(s) = K/l +-J—+Tjs) 

Although it is more tedious to implement control 
analysis for DPS using conventional techniques, DPS obey 
the same rules that a lumped system does. Those rules 

include: . . , . 

(i) The damping ratio and stability decrease with the increase 

of K p , system has a faster response; 

(ii) Integral control reduces or eliminates steady-state errors at 
the expense of reduced stability; 






(iii) Derivative control improves system stability. 

It should be pointed out that while all of these control 
techniques are familiar to most of control engineers, it is the 
structural diagram method that makes it possible to apply 
this body of knowledge to DPS. 

The step responses for proportional control where (1 ) 
k=J and two cases ofPED control where (2) K p =.01. T^l;6. 
Td-2> (3)Kp=.05, Ti=l/]0, Tj=l are shown in Fig. 3. The 
simulations were performed on the transfer functions using 
Euler's method to replace s (s-(z-l)/T) and e J =r\ where 
n=l/T. The time step used is .001 sec. 

The other classical approaches, such as lead and lag 
compensators add little significant information. 

Another more complex control problem is to regulate 
the temperature distribution along the whole x direction. 
For this case, the chosen error signal has the followinz 
form: 


e(s) = 


f 


(Q(yj) - Q*(yj))dy 


The control signal is then U(s) = W(s)e(s) , where \V(s) is 
the lumped compensator. Using the same techniques as 
above, the system response can be regulated to the desired 
specifications for a reasonably given expected signal Q*(y,s) 
(obviously, (2*(y.Jl=constant is not a reasonable choice 
because of the uniform distribution assumption of control 
signal u(t)). Fig.4 shows the steady states of the 
temperature distribution of the material being "cooled" for 
two cases of control design. 

This system will be considered further in a later paper 
where robust control methods (4] and functional observers 
will be considered. 


IV. Summary 

It has been shown that the structural theory for DPS 
has many properties similar to those of lumped systems, as 
well as some substantial differences. These are mostly related 
to having non-rauonal transfer functions. Unfortunately, 
each particular type of PDE has a corresponding non-rational 
function of the complex variable s (eg. J(xj), sinh(xj), etc). 
Thus each particular system requires unique special study to 
determine its closed loop poles as root locus methods will 
not generally be available. 
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Fig. 1 Closed-loop Poles of Proportional Control 
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Fig. 2 Closed-loop Poles of PID Control 
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Fig. 3 Step Response of Closed-loop System 



Fig. 4 Steady States of Temperature Distribution 



